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Abstract. A notion of band limited functions is considered in the case of the 
hyperbolic plane in its Poincare upper half-plane H realization. The concept of 
band-limitedness is based on the existence of the Helgason-Fourier transform 
on BL An iterative algorithm is presented, which allows to reconstruct band- 
limited functions from some countable sets of their values. It is shown that 
for sufficiently dense metric lattices a geometric rate of convergence can be 
guaranteed as long as the sampling density is high enough compared to the 
band-width of the sampled function. 



1. Introduction 

The main goal of the present article is to consider an iterative algorithm for 
reconstruction of band limited functions on the two-dimensional hyperbolic plane 
in its Poincare upper half-plane realization. 

The notion of band limited functions plays a central role in the classical signal 
analysis in which signals propagate in Euclidean space. It seems interesting to 
extend this theory to other geometries in particular to hyperbolic spaces. In this 
connection we would like to mention very interesting investigations of A. Kempf 
[10j . who used our sampling theory on manifolds to develop an approach to 
quantization of space-time and information theory. 

Analysis on the hyperbolic plane is also useful for other applied problems. As 
it was discovered by C. Bcrcnstcin and E. Casadio Tarabusi in [2], analysis on the 
hyperbolic plane plays an important role for electrical impedance imaging. An- 
other interesting applications of the analysis on the hyperbolic plane to microwave 
technology were given by A. Terras in |25j . 

The classical sampling theorem says that if / G L 2 (R) is from the Paley- Wiener 
space PW U , uj > 0, i.e., if its Fourier transform / has support supp(/) in [— 
then / is completely determined by its values at points nfi, where Q — n/uj, i.e., 



/(*) = £/(««) 



sin(w(t — nf2)) 
ui{t — nQ) 



where the equality holds in the L2-sense. The Paley- Wiener theorem states that 
/ € PWu if and only if / is an entire function of exponential type uj. 

The irregular sampling theory which was started by Paley and Wiener |14| was 
further developed by A. Beurling, P. Malliavin [3], [4], and H. Landau [13]. 
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In the early 1990's, H. Feichtinger and K. Grochenig have introduced iterative 
reconstruction methods (cf., the series of papers [5], [6], [7]) which allow to recover 
band-limited functions on R d from their irregular samples. 

On the other hand, a version of an irregular sampling theory in which band- 
limited functions are reconstructed as limits of variational splines was developed by 
I. Pesenson in the case of L2(M. d ) [20], in the case of stratified Lie groups [18], [19] . 
in the case of manifolds [17], [21], [23], [24], and in the case of a general Hilbert 
space [22j . Note that the notion of band-limited functions in Hilbert spaces was 
developed by I. Pesenson in |12j . [15] and [16] in connection with approximation 
theory and Bcsov spaces on manifolds. 

In our recent paper [8], we developed an iterative algorithm for reconstruction 
of band limited functions on general Riemannian manifolds of bounded geometry. 
In the same paper we also proposed to specify our results for manifolds which are 
of interest for applications. 

The present paper is organized as follows. In the second section we will briefly in- 
troduce Harmonic Analysis on the Poincare upper half plane H = SL(2, R) / SO(2). 
In the third section we introduce our basic notion of band limited function by using 
a non-commutative Fourier analysis on H which is given by the Helgason- Fourier 
transform. In the last section we present our iterative reconstruction algorithm. 

It is important to note that, due to the existence of rich harmonic analysis 
on H, our description of band limited functions and our reconstruction method 
become much more constructive when compared to the case of a general Riemannian 
manifold. 

2. The Poincare upper half-plane model of the hyperbolic plane 

Let G = SL(2,M.) be the special linear group of all 2 x 2 real matrices with 
determinant 1 and let K — SO(2) be the group of all rotations of R 2 . The factor 
H = G/K is known as the 2-dimensional hyperbolic space and can be described in 
many different ways. In the present paper we consider the realization of H which 
is called Poincare upper half-plane (see [9], [25]). 

As a Riemannian manifold, H is identified with the regular upper half-plane of 
the complex plane 

M = {x + iy\x,y G R, y > 0} 
with a new Riemannian metric 

ds 2 = y- 2 {dx 2 + dy 2 ) 
and corresponding Riemannian measure 

dfi = y~ 2 dxdy. 

If we define the action of a G G on z £ H as a fractional linear transformation 

a ■ z = (az + b)/(cz + d), 

then the metric ds 2 and the measure d/i are invariant under the action of G on ML 
The point i = y/—l G H is invariant for all a € K. The Haar measure dg on G can 
be normalized in a way that the following important formula holds true 

(2.1) f f(z)y- 2 dxdy= f f(g-i)dg. 
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In the corresponding space of square integrablc functions L2 (G) with the inner 
product 

< f,h>= / fhy~ 2 dxdy, 
Ju 

we consider the Laplace-Beltrami operator 

A = y 2 (d 2 x + d 2 ) 

of the metric ds 2 . 

It is known that as an operator in ^(H) = ^(H, c?/i) which is initially defined 
on Cq°(M.), A has a self-adjoint closure in L 2 (M). 
Moreover, if / and A/ belong to £ 2 (H), then 

<A/,/><-i||/|| 2 , 

where ||/|| denotes the L 2 (IHI) norm of /. 

We can introduce the full scale of Sobolev spaces H a = H a (M),a e R. For a > 
the space H a is understood as the domain of the operator (— A)"'/ 2 in the sense of 
the general theory of positive definite sclf-adjoint operators with the graph norm 

||/IU = 11/11 + ||A-/ 2 /||. 
For a < 0, H a is understood as a dual space to H~ a with respect to the paring 

<f,h>= / fhdpi. 
Ju 

The Helgason-Fourier transform of / is defined by the formula 

/(,,,)= //wiSRJFi,-**, 

Ju 

for s G C, ip € (0, 2tt], and where k v <E SO(2) is the rotation of R 2 by angle ip. 

We have the following inversion formula for all functions / from the space (H) 
of infinitely differentiable functions on EI with compact support: 

f(z) = (8TT 2 )- 1 [ [ f(it + 1/2, ip)Im(k v z) it+1/2 ttanhTTtdipdt. 

The Plancherel Theorem states that a map / — > / can be extended to an isometry 
of L 2 (H) (with respect to invariant measure d[i) onto L 2 (1R x (0, 2ir]) (with respect 
to the measure (87r 2 ) _1 t ta,nh Trtdtdip). 

If / is a function on H and ip is a K — S l O(2)-invariant function on H their 
convolution is defined by the formula 

/ * V{9 -i)= /O^ 1 • i)p{u)du, i = V^T, 

JSL(2,R) 

where du is the Haar measure on SL(2,R). It is known, that for the Helgason- 
Fourier transform the following formula holds true 



(2.2) 



f*P = f-<P- 
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3. Band limited signals and metric lattices 

In this section we introduce notions of band limited signals, metric lattices, 
develop some of their properties and prove inequalities which are used later. 

Definition 1. We will say that / £ L,2(M) is w-band limited function if f(s, tp) = 0, 
for \s\ > uj. The set of all w-band limited functions will be denoted as _E W (H). 

According to this variant of Harmonic Analysis on H, the following formula 
holds: 

(3-1) Af = — (V + 1) /. 

This formula, the Plancherel theorem, and the Sobolev embedding theorem imme- 
diately imply the next result. 

Theorem 3.1. For any real a and every w > the set E^iM) is a closed invariant 
linear subspace of the Sobolev space J3°"(M). Iff £ E U (H.) then, for every a > 0, 
the function A CT / belongs to C°°(M) and is bounded on H. 

Moreover, the same formula (3.1) and the Plancherel theorem give us the follow- 
ing inequality (3.2), which is appropriate to call the Bernstein inequality. 

Theorem 3.2. A function f belongs to the space E u (M),u> > 0, if and only if for 
every del, the following inequality holds true 

(3-2) ||AV||< (w 2 + ^ ll/ll- 

Proof. By using the Plancherel formula and (3.1) we obtain that for every uj- band 
limited function 

l|Aff/l|2 ^i S t<u jJ ( t2 + i)V(^)| 2 "anh^dt< + \ 



>\t\<u JO 

ti: 
1 



Conversely, if / satisfies (3.2), then for any e > and any a > we have 

r 2ir 

\f{t, ip)\ 2 tta,-ah-Ktd<pdt < 



8tt 



t\<ui+e JO 



(3-3) 



2,1 \ 2(T 

^ + \ 



It means that, for any e > 0, the function f(t, ip) is zero on (— uj — e, uj + e) x (0, 2tt]. 
The statement is proved. □ 

We take the point iei and consider a small open ball B(i, r/4) in the invariant 
metric ds 2 = y~ 2 (dx 2 + dy 2 ). It is possible to find such elements gj £ G that the 
family of balls B(xj,r/4) 7 Xj — gj ■ i, has the following maximal property: there 
is no ball in H of radius r/4 which would have empty intersection with every ball 
from this family. Then the balls of double radius B(xj,r/2) would form a cover of 
H. Of course, the balls B(xj,r) will also form a cover of H. Let us estimate the 
multiplicity of this cover. 
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Note that the Riemannian volume B(p) of a ball in H is independent of its center 
and is given by the formula 

rp 



B(p) = 2tt f sinhidi. 
Jo 



Every ball from the family {B{x vi r)} that has non-empty intersection with a 
particular ball B(x fl7 r) is contained in the ball B(x^,3r). Since any two balls from 
the family {_B(x„, r/4)} are disjoint, it gives the following estimate for the index of 
multiplicity of the cover {B(x Vl r)}\ 



B{ir) J Q r sinh tdt |oi2 : ,„ 
B{r/A) ~ j^smhtdt 



(3-4) ^Bk < °/ 4 < (12)^ 3 



In fact, the last estimate can be improved. 
We have proved the following Lemma. 

Lemma 3.3. For any sufficiently small r > there exists a set of points {xj} C H, 
such that: 

1) open balls B(xj,r/4) are disjoint; 

2) open balls B(xj 7 r/2) form a cover o/H; 

3) multiplicity of the cover by open balls B(xj,r) is not greater N — N(M) = 
12 2 e 3 . 



Definition 2. We will use notation Z({xj},r, N) for any set of points {xj} e H 
which satisfies the properties 1)- 3) from the Lemma 3.3 with a positive small r 
and we will call such set a metric (r, 7V)-lattice of H. 

Let us emphasize that an (r, A^)-lattice assumes a "uniform" distribution of 
points on H with respect to the hyperbolic metric. It will not look as a "uniform" 
distribution from the point of view of the Euclidean metric: it will become denser 
and denser when approaching the real line and will become sparser and sparser 
when going to infinity. 

Sobolev spaces H a (H) , a € R, can also be described in terms of local geodesic 
coordinates on H. 

In what follows we fix a small positive r > and an (r, 7V)-lattice in EI which will 
be denoted by Z({y u } 7 r, N). For the corresponding set of balls {B(y u , r)}, which 
satisfy Lemma 3.3, we consider a uniformly bounded partition of unity ip = {Vv} 
associated with {B{y Vl r/2)}. Namely, we give the following Definition. 

Definition 3. A uniformly bounded partition of unity tp = associated with 

{B{y Vl r/2)} is a set of non-negative C°°(H) functions such that 

a) suppVv C B{y v ,r/2), 

b) sup^ |<9' Q '?/v| < C(a), where C(a) is independent of v. 
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An equivalent norm on H a can be introduced by the formula 

(3.5) \\n^ (j>2w»f\\ 2 H « m ) . 

where ||/||jj°-(k2) denotes the regular Sobolev norm on the plane M 2 . 

The following important result is an adaptation of the Lemma 3.3 from |21j for 
the case when the dimension of the manifold is d = 2. 

Theorem 3.4. For any k > 1 there exist a constant C = C{k, N) > such that for 
any sufficiently small r > and any (r, N) -lattice Z = Z({x^},r, N) the following 
inequality holds true 



(3.6) H/ll < C { r [ \f(xj)\ 2 ) + ^11 A fe/2 /ll }, k > 1. 

4. Iterative reconstruction algorithm 

In this section we describe our first iterative algorithms for reconstruction of 
band-limited functions on manifolds. The following Lemma 4.1 presents a generic 
idea which is used for reconstruction of different classes of analytic functions. This 
result basically says that if a bounded operator is close to identity operator then it 
can be inverted by a Neumann series. 

The goal of section is to realize this idea in the case of band-limited functions 
on manifolds. 

Lemma 4.1. Let A be a bounded operator in a Hilbert space H such that for some 
7 < 1 and for all f e H 

(4-1) ||/_A/||< 7 ||/||. 

Then A is invertible and f can be recovered from Af by the following iterative 
procedure. If /q = Af and 

(4.2) /„ +1 = /„ + A(f - /») 
then 

(4.3) lim /„ = / 

n— >oo 

with the error estimate 

(4-4) ||/-/„||<7 n+1 ||/ll- 

Proof. Since the norm of the operator / — A is less than 1 , operator A is invertible 
and, thus, its inverse A -1 can be represented by a Neumann series A -1 = ^2 n=0 (I— 
A) n . Thus, every / G H is determined by Af, because 

oo / n \ 

(4.5) f = A- 1 AfY J (I~A) k Af lim £(J-A)Mf . 

n=0 \n=0 / 
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The sequence /„ of partial sums, /„ : X)fe=o(^~ A) k Af, satisfies the stated recursion 
relation. Indeed, 



n+l n+1 

f n+1 J2(I - *) k Af = Af + J2(I- A)Af 

k=0 fc=l 



Af + (I - A) J2(I - A) k Af = Af + A)f n = /„ + A(f - /„). 

k=0 

Moreover, the identity 

oo 

{I~A) k = (I-A) n+1 A-\ 

k=n+l 

implies that for all / G E U1 (M) we have 



\\f-fn\\ = \\ E (/-^/||||(/-^) n+1 ^- 1 ^/|| < 7" 
k=n+l 



+ 1| 



□ 



We will apply the above Lemma in the following situation. 

Recall that in the previous section in the Definition 3 we fixed a uniformly 
bounded partition of unity {Vv} subordinated to a family of balls {£?(?/„, r/2)} and 
introduced Sobolev norms by the formula (3.5). 

Consider an (e, 7V)-lattice Z({xj} 7 e, N), e < r/2, (see Definition 2). For an 
open cover by balls {B(xj, e/2)} satisfying Definition 2, we construct a uniformly 
bounded partition of unity {6j} subordinated to {B(x } ;, e/2)}, with the properties: 

a) supp^- C B(xj,e/2), 

b) sup^ Ic?'"'^! < C(a), where C(a) is independent of j. 

Note that the assumptions that e < r/2 and that the multiplicity of the cover 
{B{y Vl r)}, which was used in the definition of the Sobolev norm (3.5), is not greater 
than N = N(M), imply that each ball B(xj,e/2) has non-empty intersections with 
no more than N(M) balls of the family {B{y Vl r/2)}. 

Given a function / e H k (M), k > I, we consider an operator 

(4.6) V z ,e(f) = Ylf( x oW j ,e j G C?{B(xj,e/2)). 

3 

It will be shown that for any lattice Z = Z({x u }, e, N), e < r/2, and correspond- 
ing uniformly bounded partition of unity 9 — {9j} subordinated to the cover 
{B(xj,e/2)}, the function Vz,e(f) belongs to L 2 (B) as long as / belongs to H k (M), 
where k > 1. 

To construct the operator A we will need an orthogonal projection from L 2 (B) 
on the space of w-band limited functions E U1 (M). 

We are going to describe this projection in terms of the Helgason-Fourier trans- 
form. In order to do so, we have to introduce an analog of the classical sine 
function, 

sinirt 

sinc(t) = . 

Trt 
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Definition 4. We definee the hyperbolic sinch CJ (g), g £ SX(2,R), to be such a 
50(2)-biinvariant function on SL(2,R) whose Helgason-Fourier transform 

sinchuis) 

is 1 for \s\ < uj, and for \s\ > ui. 

Now, we introduce the operator P^ by the formula 

(4.7) P u f = f*sinch u ,feL 2 (W). 
It is clear that the operator 

P u : L 2 (H) -> E U (B) 

is the orthogonal projection on the subspace of w-band limited functions. 
The operator A will be defined as 

(4.8) A z ,ef = P u V z , g (f), 

where Vz e is defined in (4.6) and P u is the orthogonal projection from L 2 {M) onto 
E U (B). 

The next goal is to provide a uniform estimate for the norms of the operator 
I — A Z fi on the subspace E U {S) for all (e, iV)-lattices with e < r/2. 
The following lemma provides a step in this direction. 

Lemma 4.2. For any lattice Z = Z({xj},e,N),e < r/2, and any uniformly 
bounded partition of unity 9 = {8j} subordinated to B(xj,e/2), the map V z ,e{f) 
is a continuous operator from H k (M),k > 1, into L 2 (M). In other words, there 
exists a constant C = C(M, k), such that 

(4-9) \\VzAf)h 2 (M) < C\\f\\ HHm) , 

for all f e H k (U). 

Proof. According to (3.5) we have 

and 

\\^V z ,ef\\ 2 L2{B{y „, r)) < ll^£/(^ll| 2(iW )) < ^El/(^)l 2 ' 

3 3 

where supp 9j is in B(xj,e/2) and C depends on the multiplicity N = N(M). In 
other words, 

(4.io) \\VzAf)h a <m<CeV 2 r^Mx^ 

On the other hand, by a known inequality for R 2 , 

\f(y)\<C J2 e m - 1 \\f\\ Hm{B{X] , £)) ,k>l,C = C(k), 

0<m<k 

where y € B(xj,e/2), f e C°°{B(x } ;, e)), we obtain 

\f(xj)\<C sup |/(z)| < Ce-i/ll^^.,,)), 

xeB{xj,s/2) 

where k > 1, C = C(k), and then 
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v 1/2 , x 1/2 

Ei/(^)i 2 j <ce- i r£\\f\\ 2 Hk{BiXi , e)) 

The last inequality and the inequality (4.10) give that there exists a constant C 
which depends on smoothness k and multiplicity iV(H), for which (4.9) holds true. 

□ 

We will need the following Lemma 4.3 from |S]. 

Lemma 4.3. For any k > 1 there exist a constant C = C(H, k) > 0, such that for 
any small < e < r/2 and any Z({xj},e, N) the following inequality holds true 

(4-11) £H/-/^)IU(*,,e/ 2 ) < C £ ||(/ + A) fc /2 / ||,fc> 1. 

3 

Theorem 4.4. For a given k > 1 there exist a constant C = C(H, k) > such that 
for any lattice Z({xj}, e, N), with sufficiently small e > and for any f £ E^(M.) 

(4.12) ll/-^z,fl/ll < Ce(l+u 2 )^ 2 ||/||. 

Consequently, for fixed k > 1 and lo > it is true that for any e > satisfying 



e < 

one has 



(c(l+w 2 ) fc/2 ) 



11/ - A Zi9 /|| < 7II/II wfcere 7 = Ce(l+ W 2 ) fc / 2 < 1. 
Proof. For any / e E^M) we have 

ll/-^/ll = l|p w /-p w y z , e (/)ll < ||/-Vz, fl (/)||. 

Next, since the following identity holds true 

/(*) - £ *i (*)/(*)> 

3=31 

we obtain 

ii/- vx(/)ii < ii/-E^-)ii < iiE^-E^oii ^ 

3 3 3 

£n/-/fo)iu a (-B(x 3 -,e/2)) 

3 

where the sum can be estimated by using our inequality (4.11). 
We have 

||/-A z>e /||<C £ ||(7 + A) fe / 2 /||, 

for any k > 1. 

Because, for / £ E U (M), the Bernstein inequality 

||A fe /||< w 2fc ||/|| 

holds true, it yields the inequality 

\\f-A z , f\\ <Ce(l+ W 2 ) fe / 2 ||/||. 
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For a fixed u> > and k > 1, because C = C(H, k) depends only on H and k, we 
can choose 

e< (c(l+c 2 ) fe/2 )~\ 

such that corresponding 

1 = Ce{l+uo 2 ) k ' 2 

is less than 1. Theorem is proved. □ 

Combining this result with Lemma 4.1 we obtain 
Theorem 4.5. For a given k > 1 and tu > choose an e > suc/i i/iai 

e < (cr^i+w 2 )-^ 2 , 

where the constant C — C(H, fc) taken from Theorem 4-4- Then, for any (e,N)- 
lattice Z({xj},e, N) and for any corresponding operator A — Az.e in the space 

(H) define inductively, starting from 

fo = Af: 

fn+l=fn+A(f-f n ). 

Then in the space L 2 (H) the following convergence holds true 

lim f n = f in X 2 (H), 
n— >oo 

with the error estimate 

and \\f-f n \\ < 7 " +1 ||/ll, 
where 7 = Ce(l + uj 2 ) k / 2 < 1. 
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